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Abstract 



We develop an extension of Bohmian mechanics by defining Bohm-like tra- 
jectories for quantum particles in a curved background space-time containing a 
spacelike singularity. As an example of such a metric we use the Schwarzschild 
metric, which contains two spacelike singularities, one in the past and one in the 
future. Since the particle world lines are everywhere timelike or lightlike, particles 
can be annihilated but not created at a future spacelike singularity, and created 
but not annihilated at a past spacelike singularity. It is argued that in the pres- 
ence of future (past) spacelike singularities, there is a unique natural Bohm-like 
evolution law directed to the future (past). This law differs from the one in non- 
singular space-times mainly in two ways: it involves Fock space since the particle 
number is not conserved, and the wave function is replaced by a density matrix. 
In particular, we determine the evolution equation for the density matrix, a pure- 
to-mixed evolution equation of a quasi-Lindblad form. We have to leave open 
whether a curvature cut-off needs to be introduced for this equation to be well 
defined. 

PACS numbers: 04.20.Dw; 03.65.Ta; 04.62. +v. Key words: quantum theory in 
curved background space-time; Schwarzschild space-time geometry; spacelike sin- 
gularity; Bohmian trajectories; particle creation and annihilation; pure-to-mixed 
evolution. 
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1 Introduction 

This paper, part two of a two-part series on Bohmian mechanics at space-time singu- 
larities, can be read independently of part one [31]. We consider quantum mechanics in 
a relativistic space-time with fixed background metric containing spacelike singularities 
from the perspective of Bohmian mechanics (also known as pilot-wave theory), a precise 
version of quantum mechanics in which particles have world lines. I argue that, in this 
setting, a unitary time evolution for the wave function is no longer possible, and must 
be replaced by an equation for a density matrix, see ([T]) below, an integro-differential 
equation of a quasi-Lindblad form that evolves pure states into mixed states. To my 
knowledge, this equation is novel; but it is in line with an earlier proposal of Hawking 
PZ3IIE], grounded on black hole evaporation, to the effect that the fundamental physical 
evolution law should transform pure states into mixed states. 

The role of the density matrix here is unusual: Usually, density matrices represent 
statistical mixtures, or, in the case of a reduced density matrix obtained by a partial 
trace, the state description of a system that is entangled with its environment. Here, 
in contrast, the density matrix does not represent incomplete knowledge but rather 
determines the motion of the particles, a role normally played in Bohmian mechanics 
by the wave function. Still, the evolution involves information loss, as different density 
matrices at one coordinate time may evolve into the same density matrix at a later time. 

According to the singularity theorems of general relativity [191 chap. 8], a black hole 
arising from a gravitational collapse contains a singularity, which is generally believed to 
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be spacelike. As a concrete example of our general scheme we consider N non-interacting 
spin-i particles in a Schwarzschild space-time, which contains two spacelike singulari- 
ties. The terminology and notation we use is that of quantum mechanics, rather than 
quantum field theory. The evolution equation for the density matrix then contains the 
Dirac Hamiltonian and a term connected to the singular boundary of the configuration 
space, while the configuration space arises from a spacelike hypersurface bordering on 
the singularity. Since the model assumes that the space-time metric is given, it does 
not include any back reaction such as growth of the horizon or the singularity after 
swallowing particles. 

A crucial fact for the development of our Bohm-type model is that the particle 
world lines are everywhere timelike or lightlike and thus can begin but not end on a 
past spacelike singularity (hereafter, past singularity) and end but not begin on a fu- 
ture spacelike singularity (hereafter, future singularity). That is why the discussion of 
spacelike singularities is very different from that of timelike singularities. In the absence 
of other mechanisms of particle creation and annihilation, the number of particles can 
only decrease (increase) in the presence of a future (past) singularity. Since the particle 
number is not constant, we are forced to use Fock space, usually used only in quantum 
field theory but not in quantum mechanics. The natural laws for the Bohmian particles 
specify the particles' velocities and assert that any particle that hits a future singularity 
disappears while the other particles continue moving along the appropriate Bohm-like 
trajectories. Instead of merely adding Bohmian trajectories to known "orthodox" quan- 
tum theories, we use Bohmian mechanics to find the appropriate evolution equation for 
the density matrix. 

The quasi-Lindblad equation that we obtain for the density matrix arises also in a 
different context, replacing the singularity 5? by a spacelike hypersurface S?: it arises 
from a unitarily evolving wave function if) t by tracing out those degrees of freedom 
localized in the future of . In this scenario, the density matrix does represent only 
partial information about the true quantum state, and the quasi-Lindblad equation 
represents the procedure of continuously tracing out (continuously in time) more and 
more degrees of freedom (corresponding to larger and larger portions of space). 

It might seem that there is the following alternative to our quasi-Lindblad equation: 
We may refuse to replace the wave function ip t with a density matrix, stick to the Dirac 
equation for ip t , and just accept that the time evolution is not unitary. This amounts 
to deleting the amount of wave function that has crossed the singular boundary of 
configuration space. The fact that will shrink with increasing t may seem natural 
as \\ipt\\ 2 represents the probability that no particle has hit the singularity up to time t. 
But what that really means becomes clear from the Bohmian point of view: it means 
that when one particle hits the singularity, all particles disappear. And that is a much 
less natural dynamics than postulating that when one particle hits the singularity, all 
other particles continue moving along Bohm-type trajectories. 

Due to limits of my knowledge, I have to leave open the question whether a curva- 
ture cut-off needs to be introduced to make the quasi-Lindblad equation well defined; 
I conjecture that no such cut-off is necessary. Such a cut-off can be implemented by 
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choosing an arbitrarily thin neighborhood of the singularity whose surface is a space- 
like hypersurface J^, and let 5? play the role of 5? in the quasi- Lindblad equation. It 
seems that this question boils down to the question whether the probability current as- 
sociated with a solution of the 1-particle Dirac equation possesses a continuation on the 
singular boundary I conjecture that the answer is yes. Furthermore, for the precise 
formulation of the quasi-Lindblad equation it is relevant whether the spin spaces can be 
defined on the singular boundary 5? in such a way that also a wave function obeying 
the 1-particle Dirac equation possesses a continuation on 5? . In this paper, however, 
we will simply assume that limits on 5? exist whenever needed. 

Further questions arise from past singularities as in the white hole of the Schwarzschild 
space-time. By time reversal symmetry, the same evolution equation for the density ma- 
trix that holds at a future singularity should apply here backwards in time. When looked 
at in the ordinary time direction, from the past to the future, then the evolution of the 
density matrix will be indeterministic because of the "information loss" property in the 
other time direction. This kind of indeterminism is distinct from the quantum inde- 
terminism (as expressed by the Heisenberg uncertainty relation and indeed compatible 
with the determinism of Bohmian mechanics) and the indeterminism represented by the 
stochastic law governing the configuration in Bell-type quantum field theories [HI ITT] . 
Once the evolution of the density matrix is fixed, the evolution law for the configura- 
tion is fixed by time reversal symmetry, and turns out to involve, besides a Bohm-type 
equation of motion, a stochastic law for particle creation at the singularity. 

Another work on Bohmian mechanics and black holes is Valentini's [32], proposing 
that the equivariance of the |"0| 2 distribution might be violated in the presence of black 
holes. 

This paper is organized as follows. In the remainder of this section we give an 
overview. In Section [2] we recall Bohm's law of motion for the Dirac equation. In 
Section [3] we recall the Schwarzschild metric. In Section H] we develop our new version 
of Bohmian mechanics in the presence of a future singularity, and in Section [5] in the 
presence of a past singularity. 

1.1 Horizons 

Let me make a few remarks about the status of horizons in Bohmian mechanics. 

Consider a many-particle quantum system in a background space-time containing 
a black hole. From the point of view of orthodox quantum mechanics, it is natural to 
trace out all degrees of freedom inside the event horizon. From the Bohmian viewpoint, 
in contrast, this is not natural. Instead, it is natural to ask what happens behind the 
horizon. This is so because from the orthodox viewpoint, the most important question is 
what an observer will see, and it is a frequent assumption that the relevant observers sit 
at infinity. From the Bohmian viewpoint, the most important question is what actually 
happens. Thus, to define Bohmian mechanics in a curved space-time, we need to define 
as well the trajectories inside the black hole. (A viewpoint that, like the "Copenhagen" 
view of quantum mechanics, dismisses any theory of particle positions because it regards 
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the latter as "hidden variables," may naturally tend to dismiss as unreal also everything 
that is hidden behind a horizon.) 

Thus, the Bohmian viewpoint leads us to the following attitude: What happens 
inside a black hole can and should be described by a physical theory. 

Indeed, by the nonlocal character of Bohmian mechanics, the velocities of the parti- 
cles outside the black hole may depend on the positions of the particles inside the black 
hole. But all this requires no additional research for the definition of the theory, as the 
Bohm-type law of motion we use (see Eq. (121 j) below) automatically implies influences 
across event horizons. What requires further work is not the presence of a horizon, but 
the presence of a singularity. (The need for this further work would disappear if none of 
the hypersurfaces belonging to the time foliation bordered on the singularity. However, 
for our purposes that is the uninteresting case. Furthermore, I see no reason why the 
time foliation should avoid the singularity.) 

Now we can appreciate the differences between the approach of this paper and that 
of Hawking [17] : Hawking traces out what has passed the horizon, while we trace out 
what has hit the singularity; for Hawking, taking the partial trace is only a matter of 
representing the knowledge of observers at infinity, while for us it is part of defining 
the true particle trajectories; Hawking uses positivist arguments, while we start from a 
realist picture; for Hawking, a pure-to-mixed evolution fundamentally occurs only when 
the black hole evaporates, while for us it occurs continuously during the existence of 
the singularity; for Hawking, the pure-to- mixed evolution occurs because late (i.e., post- 
evaporation) observers cannot measure early (i.e., pre-evaporation) observables, while 
for us even a demon knowing all facts has to apply the pure-to-mixed evolution; Hawking 
focuses on the quantum state at infinity (as in scattering theory), while we need the full 
time evolution. 

1.2 Background 

Bohmian mechanics was developed as a realist version of nonrelativistic quantum me- 
chanics [6] and succeeds in explaining all phenomena of quantum mechanics in terms of 
an objective, observer-independent reality consisting of point particles moving in space; 
see [2] for an overview. Bohmian mechanics possesses a natural extension to relativistic 
space-time if a preferred foliation of space-time into spacelike hypersurfaces (called the 
time foliation J 7 ) is granted [8]. This extension has also been formulated for curved 
space-time geometries (29J, 150] . but not yet for geometries with singularities. While 
horizons present no difficulty, singularities require further work to define the theory: 
Basically, we have to specify what happens when a particle hits the singularity, since 
at this point the law of motion is no longer defined. The possibility we consider here 
is that the particle gets annihilated: that is, if the system consisted of N particles, 
then the world line of the particle that has arrived at the singularity ends there, while 
the other N — 1 particles, which are not at the singularity and thus have no reason to 
vanish, continue to move according to Bohm's law of motion. To make this possible, we 
need wave functions from Fock space, i.e., superpositions of different particle numbers, 
as always when particles can get created or annihilated. Further considerations then 
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naturally lead us to specific equations, defining a Bohm-type theory 

Space-time singularities are points on the boundary of space-time where the met- 
ric cannot be extended smoothly because the curvature becomes infinite. There is no 
universally accepted mathematical construction of these boundary points from a given 
metric (see chapter 8 of [19j for a discussion), but readers may adopt the construction 
of [13], called the causal boundary, that defines boundary points as terminal indecom- 
posable past sets (TIPs) or terminal indecomposable future sets (TIFs), together with 
suitable identifications. A boundary point that does not lie at infinity is considered 
a singular boundary point. The singular boundary is timelike at a point x if x has a 
non-empty TIP and a non-empty TIF; it is future spacelike if x has empty TIF and past 
spacelike if it has empty TIP. Like a hypersurface, a singularity can be timelike at some 
points and spacelike at others. 

We will not construct the singular boundary but regard it as given. More precisely, 
we assume that space-time j$ is a manifold with boundary (211 [T], where the boundary 
represents the singularity. In particular, we assume that the singularity has the struc- 
ture of a 3-dimensional manifold. In the Schwarzschild geometry, for example, ^ is 
diffeomorphic to [—1,1] xtx§ 2 ; the singular boundary is everywhere spacelike and has 
two connected components each of topology 1 x § 2 ; one is a future singular- 

ity, the other a past singularity. We shall recall some details about the Schwarzschild 
geometry in Section 

1.3 Motivation 

It is a natural part of the research program on Bohmian mechanics to extend the theory 
to more general quantum theories, to all kinds of settings. To the extent that we have 
reason to believe that singularities exist in our universe, we obtain here a more appropri- 
ate version of Bohmian mechanics. Concerning techniques of constructing Bohm-type 
models, we find that the Bohm-type law of motion proposed by Diirr et al. [8] for rel- 
ativistic space-time with a foliation works unproblematically even in the presence of a 
singularity, a result that lends support to this law of motion. 

Since Bohmian mechanics is a particularly precise and unambiguous version of quan- 
tum mechanics it may serve as a tool for studying quantum mechanics in curved space- 
time. Thus, our study can as well be regarded as one on quantum mechanics at space- 
time singularities. 

It is remarkable that the presence of spacelike singularities forces us to change the 
structure of Bohmian (or quantum) mechanics in three ways: First, as a consequence 
of the non-conservation of particle number, it is necessary to use Fock space; second, a 
density matrix replaces the wave function in its role of guiding the particles; third, the 
evolution is no longer unitary. 

Independently of the Bohmian viewpoint, we obtain a novel evolution equation for 
the density matrix in the presence of a spacelike singularity. It can be argued that 
this evolution of the density matrix is more fundamental than the unitary dynamics 
of quantum mechanics, since if our universe contains any spacelike singularities then a 
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unitary evolution apparently does not exist. Conversely, if we think that the fundamental 
evolution should be unitary, we should be skeptical about the existence of spacelike 
singularities. In contrast, timelike singularities do not enforce deviations from unitarity 
and are thus less dramatic in Bohmian (or quantum) mechanics than spacelike ones; this 
situation differs from that in classical mechanics, where future spacelike singularities are 
unproblematic while timelike and past spacelike singularities are not covered by the laws 
of classical mechanics — anything could come out of such singularities. 

When using a foliation that consists of Cauchy hypersurfaces (which, in particular, 
do not border on the singularity), the usual unitary evolution works and no need for 
a quasi-Lindblad evolution arises. So why even consider a non-Cauchy foliation? One 
reason is that for a non-Cauchy foliation, something unusual and interesting happens: 
the unitary evolution gets replaced by the quasi-Lindblad evolution. A second and more 
serious reason is that, according to Bohmian mechanics (in the version that we use here), 
there is a fact in nature about which foliation is the right one, i.e., which is the time 
foliation T . Since the law selecting T is presently not known, we have to admit the 
possibility that T is a non-Cauchy foliation. Indeed, a candidate for a law for T is 
formulated in (1241) below, and this law leads inevitably to non-Cauchy foliations in the 
presence of a singularity. This observation urges us to consider non-Cauchy foliations. 



1.4 Overview of New Equations 

Let us have a first look at relevant equations, postponing a more careful definition 
and deeper discussion to Section @J Let y be a future singularity, regarded as a 3- 
dimensional manifold, the boundary of space-time. 

The evolution equation for the density matrix p t on (either bosonic or fermionic) 
Fock space is an integro-differential equation and reads 

^ = ^p t ,H t }+£p t , (1) 

with H t the Hamiltonian on Fock space arising from the Dirac equation; the square 
brackets mean the commutator. The symbol £ denotes a particular superoperator (i.e., 
an operator acting on density matrices rather than wave functions) which can be defined 
explicitly by 

(Cp t )(q;r) = V(#g + l)(#r + 1) J d 2 xw(x) x 

x lim c^x) d^y) tr Sy (p t (q,y;r,y) a^(y) ) . (2) 

Let us explain the notation. Let E t be the spacelike hypersurface corresponding to time 
t and S y the spin space at space-time point y; with a configuration q = (qi, ... , q^) 6 
there is associated the spin space 

S q = S qi <g> • • • ® S qN . (3) 
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We write density matrices p as spinor- valued functions p(q; r) of two configuration vari- 
ables, corresponding to the formal expression 



p(q-r) = (q\p\r) , (4) 

where the scalar product is a partial scalar product taken only over the position degrees 
of freedom but not the spin degrees of freedom; p(q; r) is a linear mapping S r — > S q . The 
symbol j^q denotes the number of particles in the configuration q. Let 5^ t = H S t , 
n^ix) the unit normal vector on S ( at x G S t , and 

a^{x) = (n v {x) 1 u {x)y l r{x), (5) 

where 7 M are the same gamma matrices as in the Dirac equation. Furthermore, in ([51) 
d 2 x denotes the coordinate area measure of the surface element; w(x) dt is the thickness 
in coordinate space of the 3-dimensional strip between the 2-dimensional surfaces S^t 
and S^t+dt at x G S?t\ c n{ x ) is the 4- vector which in coordinate space is orthogonal to 5? 
at x, points to the future, and has unit length; <i 4 is the determinantal 4-volume factor, 

<k{y) = \/-detg^(y) ; (6) 

trs y denotes the partial trace over S y ] and the pair q,y means the configuration arising 
from q by adding a particle at yu 

The Bohm-type trajectories are defined by a first-order equation for the particle 
configuration Q(t) = (Qi(t), . . . , Qn^)) G S 



jv 
t ' 



dQt d A {Q k ) tT Sg (p(q;q)aZ(q)) 



dt d 3 (Q k ) tr Sq p(q;q) 
where d% is the determinantal 3-volume factor 



(7) 

q=Q(t) 



d 3 (x) = y/-detgW(x) (8) 

and 

o£( q ) = I®( k -V ® a^q k ) ® ^(N-k-i) ? (9) 

for which we will often simply write a M (gfc), as the argument q k makes it unambiguous 
on which spin index the matrix acts. 

In addition, it is postulated that whenever the configuration Q(t) reaches the bound- 
ary of configuration space, which happens when one of the particles reaches the singu- 
larity, then that particle will be annihilated, corresponding to a jump 

Q{t-)^Q{t+)=Q{t-)\y (10) 



1 Although some quantities in ([2]) explicitly refer to a coordinate system, the formula ([2]) is actually 
equally valid in every coordinate system that has the t function as its time coordinate. Concern- 
ing the limit y — > x,y ^ & ', I conjecture that the spin bundle can be extended to the singularity 
in such a way that pt(q,x;r,x) is well defined for x G 5? without the need of a limit, and that 

lim c M (x) di (y) a M (y) is a well-defined endomorphism ay (x) of S x . 
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in the configuration, where Q(t±) = lim £ \ Q(t ± e). 
The equation of motion (j7|) is of the form 

dQ(t) j(Q(t)) 



(11) 



dt ' p(Q(t)) 

with the probability density p (relative to coordinate volume) given by 

p(q) = ppt(q) = d 3 (q) tr Sq p t (q; q) , 



(12) 



where 



(l3(?) = nV-det#( 9fc ), 




(13) 



k=l 

and the probability current density j given by 



3k fa) = fk^il) = d z(l \ Qk) di{q k ) ti Sq (p(?J q)oc^{q k )) . 



(14) 



Indeed, (J7J) is the straightforward generalization of Bohm's equation of motion from 
wave functions to density matrices; by setting p = \4>)(?p\ we obtain the usual form of 
Bohm's equation of motion for the Dirac equation, 



with (ffip = 07°^ the inner product in spin space. The version (£7§ with a density matrix 
has been considered before in two places: 

• in [3] for contrasting the trajectories arising from (J7]) with those arising from 
(|15p for a random ip whose distribution has density matrix p, illustrating that in 
Bohmian mechanics nature needs to know ip, not just p; 

• in [TD] for finding a replacement of the conditional wave function for particles with 



Since, according to d7j), the density matrix p has the role of "guiding" the particles, 
to determine how they move, it is a fundamental density matrix (in the terminology of 
[10J), as opposed to a description with less than full information. Eq. defines a "pure- 
to-mixed" evolution, where a density matrix is called "pure" if it is a one-dimensional 
projection and "mixed" otherwise. Readers should not understand "mixed" as referring 
to any ensemble of wave functions, or as indicating the presence of any randomness. 
The density matrix plays not a statistical role here, but rather the role of a fundamental 
object in the theory; this difference is clear from the Bohmian perspective but hard to 
discern from the perspective of orthodox quantum theory. 

The evolution equation ([T]) is of a quasi- Lindblad form, as we explain now. The 
Lindblad equation reads [221 EH] 



dQt _ d,(Q k ) riQh^QkMQ) 



(15) 



dt d 3 (Q k ) ip*(Q)ip(Q) 



spin. 




(16) 



A 



A 
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where A x is a countable family of bounded operators and H a bounded self-adjoint op- 
erator, and the curly brackets mean the anti-commutator. Every (uniformly continuous) 
quantum dynamical semigroup is generated by an equation of this form [221 HB] • A more 
general version reads 

^ = {[p u H t ] + f iit(d\)A t (\)p t A* t (\) - \{p t ,J p t (dX)A* t (X)A t (X)} , (17) 

where A is any parameter, p t any measure over that parameter, A t (X) any operator, and 
H t a self-adjoint operator. 

Our equation ([1]) is similar but not exactly of this type. Rather, it is of the form 

^i = i[p u H t } + J p t (dX) A t (X) p t A*(X), (18) 

and that is why we call it a quasi-Lindblad equation. There are two main differences 
to (TTTj) : First, our Hamiltonian H t is not self-adjoint, a fact corresponding to the loss 
of probability at the singular boundary of configuration space. Second, the third term 
on the right hand side of (fT7|) is missing in ([TBI . To understand why, note that that 
term is needed in ( |17p to compensate the trace of the second term: if the trace of p t is 
to be conserved, the trace of the right hand side of (1171) should vanish; a commutator 
[p t , H t ] with a self-adjoint H t has zero trace; since the second and the third term contain 
the same factors cyclically permuted, they have equal trace — but opposite signs in front 
of them. In contrast, the trace of the commutator in ([TBI , or ([TJ), is not zero, but 
instead minus the rate of probability loss at the boundary of configuration space. It is 
equally large, but with opposite sign, as tr Cp t , so that also the right hand side of ([TBI 
is traceless. 

In order to show that ([1]) is of the form (TlB]) we need more details about how the 
Hilbert space is defined; we postpone this to Section H~Tl see equation (IM1) . 

For t 2 > ti, the mapping p tl \— > p t , 2 defined by solving (pQ) can be regarded as a super- 
operator S l t l from the trace class of M[ x to the trace class of The superoperators 
are completely positiv^] and satisfy S^p = p and S^S^ = S^. This structure is similar 
to that of quantum dynamical semigroups, which often arise as effective evolution of 
reduced density matrices, for example representing decoherence. Here, of course, it is 
not decoherence that causes the pure-to-mixed evolution; it is not a reduced density 
matrix that obeys ([TJ but the fundamental density matrix; and ([T]) is not an effective 
equation but the fundamental law of nature replacing the Schrodinger equation. 

The question arises how it can be possible that the equation (JlJ fails to be of the 
Lindblad form (fT6j) or (|T7|) . in view of the Lindblad theorem [22l [T6] that (uniformly 
continuous) quantum dynamical semigroups are always generated by an equation of 
that form. One might suspect the reason is that the superoperators S\l defined by 

2 A linear mapping from the trace class of 3f\ to the trace class of is called completely positive if, 
first, it maps positive operators on ,¥P\ to positive operators on M2 and, second, its obvious extension 
maps positive operators on Jfi ® j£| to positive operators on ® ^3 , for any ^ [7] . There is wide 
consensus that this is the appropriate positivity property for transformations of density matrices. 
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(PQ) do not form a semigroup, since the semigroup structure requires time-translation 
invariance (in the sense of Sq = S l s +S ). While it is true that the do not usually form 
a semigroup, one can devise cases that actually are time-translation invariant and for 
which the do form a semigroup (using a suitable identification of the Hilbert spaces 
associated with t\ and £2); so the semigroup assumption is not the relevant hypothesis 
of the Lindblad theorem that fails here. What fails is only the assumption of uniform 
continuity, so I conclude that this assumption is not, as one might have thought, merely 
a technical assumption without physical relevance, but a strong restriction excluding 
quasi-Lindblad equations. 

For past singularities, the evolution of the density matrix into the future is not 
uniquely determined. If we fix one such evolution, we obtain a unique law for the rate 
at which particles are created at the singularity: The rate of creation of a particle 
at time t in the 2-dimensional surface element d 2 x in ynSt, given that the present 
configuration is Q(t) = q, is 

w(x) (#g + 1) lim c^x) d A (y) tr Sq y (p t (q, y; q, y) a"(y)) 

= — tr^w) ^ < 19 > 

This formula is closely related to the creation rate formula of so-called Bell-type quan- 
tum field theories [91 [HI HL versions of Bohmian mechanics with particle creation and 
annihilation. 



2 Bohmian Mechanics 

Bohmian mechanics [61 HH, EH] is well understood in the realm of non-relativistic quan- 
tum mechanics, but needs further development in the directions of relativistic physics, 
quantum field theory, and quantum gravity. This paper concerns the relativistic exten- 
sion in a classical gravitational field, but connects also with quantum field theory. 

Bohmian mechanics postulates that particles have trajectories, governed by an equa- 
tion of motion of the type 

dQt j+{Qt) . , 

dt p*(Q t )' 1 ) 

where Q t is the position of the particle at time t (or, for a system of several particles, 
the configuration) , and y and p^ are, respectively, the quantum mechanical probability 
current and probability density as determined by the wave function ip. As a conse- 
quence of the structure ( 1201) of the law of motion, if at any time t the particle position 
(or configuration) is random with distribution p^*, then this is also true of any other 
time t. This property is called equivariance. As a (quite non-obvious) consequence of 
that, inhabitants of a Bohmian universe, consisting of these particles with trajectories, 
would observe the same probabilities in their experiments as predicted by the quantum 
formalism [12]. That is how Bohmian mechanics explains quantum mechanics. In fact, 
Bohmian mechanics accounts for all phenomena of non-relativistic quantum mechanics. 
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2.1 In Relativistic Space-Time 

With the invocation of a preferred foliation T of space-time into spacelike hypersurfaces, 
given by a Lorentz invariant law and called the time foliation, it is known [8j [29j [30] that 
Bohmian mechanics possesses a natural generalization to relativistic space-time. The 
possibility of a preferred foliation seems against the spirit of relativity, but certainly 
worth exploring. It is suggested by the empirical fact of quantum non-locality, and it 
is suggested by the structure of the Bohmian law of motion (I2"0"j) for many particles, 
in which the velocity of a particle depends on the instantaneous position of the other 
particles. Using a time foliation JF, a Bohm-type equation of motion was formulated in 
[8] for flat space-time, and the straightforward generalization to curved space-time was 
formulated and mathematically studied in 



dQ 



i i k 



ds 



cc (Q 1 (E), Q N (E)) J] n w (Q f (E)) , (21) 



where Qk{s) is the world line of particle k G {1, . . . , N}, s is any curve parameter, E is 
the hypersurface in T containing Qk{s), n(x) is the unit normal vector on E at x G E, 
Qi(E) is the point where the world line of particle i crosses E, and 

jui-m = ^yii g,...® yw)^, (22) 

is the probability multi-current associated with the iV-particle Dirac wave function if). 
This wave function if) could either be a multi-time wave function defined on ^ , where 
J% is space-time, or, since we never use if) for configurations that are not simultaneous, it 
suffices that if) is defined on the 3N + 1-dimensional manifold UseJ 17 ^ N °^ simultaneous 
configurations. The probability density, relative to the invariant volume on S, is given 
by 

N 

j^ N (qi,...,q N ) ]JnM (23) 
fc=i 

for q x , . . . ,q N G E. 

The extension of Bohmian mechanics to relativistic space-time that we just described 
does not automatically include, however, space-time geometries with singularities. The 
treatment of singularities requires some fundamental extensions of Bohmian mechanics, 
and forms a test case for the robustness of the equation of motion (I2TI) . Well, the 
equation has stood the test, both with timelike [3TJ and spacelike singularities. 

As mentioned before, the foliation might itself be dynamical. An example of a 
possible Lorentz invariant evolution law for the foliation is 

V^n v - V v n^ = , (24) 

which is equivalent to saying that the infinitesimal timelike distance between two nearby 
hypersurfaces from the foliation is constant along the hypersurface. This law allows to 
choose an initial spacelike hypersurface and then determines the foliation. A special 
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foliation obeying fl2^|) is the one consisting of the surfaces of constant timelike distance 
from the big bang. Note, however, that the law of motion (|2"TT) does not require any 
particular choice of law for the foliation, except that the foliation does not depend on 
the particle configuration (while it may depend on the wave function). Note further that 
in a space-time with horizons and singularities, a foliation law like (12^|) will frequently 
lead to hypersurfaces lying partly outside and partly inside the horizon, and indeed to 
hypersurfaces bordering on a singularity. 

2.2 Adapted Coordinates 

When expressed in terms of coordinates that are adapted to the time foliation, i.e., such 
that the time coordinate function is constant on every time leaf S G J 7 , ( 1211) assumes 
the form ([T5]) . while the probability density is given by 

p(q) = d 3 (q)r(q)H<l), (25) 

and the current by 

fM) = d 3 (q \ q k ) d 4 (q k ) lf;*(q) a^q k ) lfj(q) . (26) 

These equations need some elaboration. 

Let us first turn to the definition of the inner product. At x G except on the 
singularity, every future-pointing timelike vector n M G T x j& defines a positive definite 
Hermitian inner product in the spin space S x , usually denoted 

?Vy"V (27) 

for any <f>,if) G S x . We will always use the inner product in S x defined by the future- 
pointing unit normal vector n M = n^(x) on the unique time leaf E G T passing through 
x, and denote that inner product by <p*ip for any <j>, if) G S x . In terms of this inner 
product, 

= (f)* a^ip (28) 
with given by (jSJ); both 7 and a are cross-sections of the vector bundl^f] 

[J CT X JK ® End(S x ) , (29) 

where CT X ^ denotes the complexified tangent space and End(S x ) the space of endo- 
morphisms of S x . From ( J28l) it follows that 

jui-m = ^*(q,mi g) . . . ® (30) 

3 We use the notation U x( zbE x , rather than E — > B, to denote the vector bundle over the base 
manifold B with fiber spaces E x , assuming it is clear from the context which bundle structure (as 
defined by the bundle maps) and, if appropriate, connection is intended. 
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and 

n^x)a^x)=I, (31) 

where I means the identity. 

Expressed relative to an orthonormal basis in T x ^ with as the timelike basis 
vector, and relative to the associated basis in spin space S x , the gamma matrices assume 
their standard form; the basis in spin-space is orthonormal, <p*ip = 5^ s= i(0s)*V's f° r t ne 
components 4> s ,ip s G C of 0, ip G S x ] a is the identity; and a 1 , a 2 , a 3 are the standard 
Dirac alpha matrices. 

More generally, in any coordinate system adapted to J 7 , a is a multiple of the 

identity, namely 

a°(x) = y/g™(xjl. (32) 

To see this, note that for every vector u G T x ^, 

u° = dt(u) = vV (x)rv(x) , (33) 

where dt denotes the 1-form obtained by differentiating the t function. As another 
consequence of (133]). we can re- write the right hand side of f[2~TT) in adapted coordinates 
as 

•0...0,A*,0~0(Q( t )) TT 1 . (34) 

If we parameterize the world lines by the time coordinate then dQ\jdt = 1, so we can 
re-write (1211) as 

dQl f-^-°(Q(t)) 



dt f-°(Q(t)) 
with the index [i in the fc-th place. 
We also note the formula 



(35) 



which follows from the fact that d±{x), the Lorentzian 4- volume spanned by the coordi- 
nate basis of T xt y#, equals the 4- volume spanned by the future-pointing vector w G T X ^M 
that is normal (in the sense of g^) to S t and has w° — 1, together with the 3 space- 
like coordinate basis vectors. Due to orthogonality, this 4-volume is the product of the 
Lorentzian length of w and the Riemannian 3- volume spanned by the 3 spacelike basis 
vectors, which is d 3 (x). That is, d 4 (x) = ^w^w u g^ iu d 3 (x), while = (V iy tV^t)~ 1 V At t, 
which implies ( 136|) . 

As a consequence of (|32|) and (|36|) . 

d 4 (x)a°(x) = d 3 (x)I. (37) 

Now we are ready to determine the probability and current density. While the tensor 
j'w-mjv re f er s to invariant volume, we prefer to express all densities relative to coordinate 
volume because usually invariant volume, but not coordinate volume, becomes singular 
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at a space-time singularity. We thus obtain from (1231) the formula (1251) for the probability 
density. The corresponding formula f[2"B"l) for the current can be derived as follows. 
Regarding f[2"5l) as known, we have j®(q) because it must equal the density p(q); to see 
that this agrees with (1261) . use (1371) . The spacelike components can be obtained from 
the fact that the quotient of the current and the density is the velocity, and thus must 
be the same as in (1351) : 

p(Q) dt f-°(Q) • 1 ) 

As a consequence of ( 1251) and ( 1371) . 

p(q) = dMf-°(q) (39) 

with 



^(g) = n v ~ det g ^ qk) • ^ 40 ^ 

k=l 

Hence 

^(?) = ^(g)/- w - Q (41) 

with n in the fc-th place, which implies (1261) . 

3 Schwarzschild Space-Time 

As an example space-time we use the Schwarzschild space-time, which we take to 
be Kruskal's maximal extension of the Schwarzschild metric EDI ESI El], and which 
features two spacelike singularities, one in the past and one in the future. 

3.1 Definition 

We use the Kruskal coordinates t' , x', ip, in which the metric is given by 

ds 2 = F 2 dt' 2 - F 2 dx' 2 - r 2 (d$ 2 + sin 2 dip 2 ) (42) 

with 

F 2 = ^Ml e -^M (43) 
r 

and r determined implicitly by the equation 

e - x ' 2 = -{r -2M)e rl2M . (44) 
The coordinates t' and x' only take such values that 

t' 2 - x' 2 < 2M . (45) 

We note for later use that 

d 4 (t', £',$,¥?) = FVsintf, d 3 (t',x',>&,(p) = Fr 2 sin$. (46) 
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The sing ularity lies at t' 2 - x' 2 = 2M . That is, the space-time ^£ is the manifold- 
wit h-boundary given by 

J[ = {(if, x ') E R 2 : t' 2 - x' 2 < 2M} x § 2 , (47) 

and the singularity 5? = is 

5P = { (t\ x') : t' 2 - x' 2 = 2M} x § 2 , (48) 

which has two connected components, 5? — 11=5*2, with 



y x = {(£', x ') : t' = V2M + x' 2 } x § 2 (49) 

and 

y 2 = {(t', x ') : t' = -V2M + x 12 } x § 2 . (50) 
is a future singularity, and ^ 2 a past singularity 

A curvature cut-off can be implemented by cutting out from ^# a neighborhood of 
the singularity, thus making a spacelike hypersurface the new boundary. For example, 
in order to cut out from the Schwarzschild space-time we could set 



S? = {{f, x') : t' = V2M + x' 2 - £ }x§ 2 (51) 

with e>0a small constant. 

As the time foliation T we take the level surfaces of the t' function. Note that these 
hypersurfaces are not Cauchy surfaces, except for —\J2M <t'< \[2M. For t' > \[2M 1 
E t / borders on the singularity Jt^i, as D Tif ^ 0, and S t / consists of two connected 
components0 corresponding to x' > y/ 1' 2 — 2M and x' < —\/t' 2 — 2M, each of topology 
[0, oo) x § 2 . For t' > V2M, the singular boundary dYj t < consists of two unconnected 
spheres, one corresponding to x' = \Jt 12 — 2M, the other to x' = —y/t' 2 — 2M; for 
t' = \J2M the two spheres coincide. Likewise, for t' < — \/2M, E t / borders on the 
singularity and consists of two connected components, each of topology M 3 minus 
an open ball; for t' < — V2M, the singular boundary <9£ t / consists of two unconnected 
spheres. 

Among other spacelike foliations, there are some whose leaves will, like those of J 7 , 
border on the singularity and others whose leaves will not, e.g., T = {S s : — 1 < s < 1} 
with 

S s = {if = sV2M + x' 2 }. (52) 

Indeed, all of these leaves are Cauchy surfaces; as a consequence, such foliations will be 
uninteresting to us, as the time evolution is unitary and the pure-to-mixed evolution we 
introduced does not occur. 



4 Needless to say, particles in different components can be entangled with each other. As a conse- 
quence, the Bohmian velocity of one particle may depend on the position of the other, and results of 
experiments carried out in different components can be nonlocally correlated. 
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3.2 End Points of Causal Curves 



Finally, we mention that every causal curve (i.e., one that is everywhere timelike or 
lightlike), if it is future inextendible and does not reach infinity, has an end point on 
the singularity In particular, it is not possible that the curve has more than one 
accumulation point on 5^\. 

To see this, consider such a curve x(t') that cannot be extended beyond time t' ; to 
see that the x' coordinate converges as t' — > t' Q , note that for any e > and t' Q — e =: 
t[ < t' < t' , x'(t r ) must lie between x'(t[) — e and ^'(i'J + e. For showing that also 
the angular coordinates converge it suffices, by rotational symmetry, to consider the <p 
coordinate and show that the function <p(t') has bounded variation in the open time 
interval (i'i,to). The total variation of <p(t') in this interval is 



V 



rt'o 

/ dt' 


dip 


It'i 


dt' 



(53) 



Since the curve is causal, ds > 0, we have from fl42ft that 



dip 



dt 1 



F 
< — 

r 



AM 



-r/4M 



(54) 



We find it useful to parameterize the curve by r rather than t'; to this end, we obtain 
from (jUJ (noting that dr < 0) that 



dt' 



re 



r/2M 



2M(2t' + 2x> dx'/dt 1 



-dr < 



re 



r/2M 



4M(t[ - \x'(t[)\ 



-dr 



(55) 



using \dx'/dt'\ < 1 and choosing e so small that t! + x' dx'/dt' > t[ — — s > 0. 

Inserting (1341) and (1531) into (|53|) and using r > 0, we obtain that 



V < 



< 



HA) 4M 
e r(t[)/4M 



-r/AM _ 



re 



r/2M 



(t[-\x'(t[)\ 
which completes the proof. 



4M(t; - \x'(t[)\ -e) 

r(t[) dr e r(t[)/4M 



r (^-IxTOI-e) 



2^r(t[) < oo 



(56) 



4 Future Spacelike Singularities 

4.1 Mathematical Framework 

We now set up the mathematical structure; since we do not strive for mathematical 
precision, not every concept will be sharply defined. 

We assume that we are given a 4-manifold with boundary ^# as the space-time 
(where 5? = cL# is the boundary), equipped with a Lorentzian metric g^ u on the interior 
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^tft° = ^ \ 5? , a time orientation, a complex vector bundle over ^ of spin spaces which 
we denote by S = U x£ ^S x , and a foliation into spacelike hyper surf aces, the time foliation 
T . We assume that 5? = <9^# is a future singularity. For simplicity, we assume further 
that T is parameterized by a time parameter t, i.e., T = {E t : t\ < t < t^}, and that 
this is done in such a way that t as a function on ^# (defined by t(x) = r <^ x G E T ) 
has nowhere-vanishing gradient, Vx G ^ : Vt(x) 7^ 0. For the Dirac equation, S has 
complex rank 4 over 
Let 

00 

Q = Q(E) = |J E" (57) 

?1=0 

be the configuration space of a variable number of particles associated with the hy- 
persurface E; the union is understood as a disjoint union; elements of Q are ordered 
configurations, but the ordering is physically irrelevant. We write j^q for the number 
of particles in the configuration q. On the set Q we consider two measures, the Rie- 
mannian (invariant) volume d 3 (q) dq and the coordinate volume dq, with d 3 as in fll3p . 
They arise from two measures on E, the Riemannian (invariant) volume d 3 (x) dx and 
the coordinate volume dx, with d 3 (x) = y — det (x). The measure d 3 (x) dx is the 
one that arises from the Riemannian metric on E inherited from g uu . 

The manifold Q is further equipped with the vector bundle U q£ QS q defined in (j3J) from 
the bundle U x( zy,S x . Throughout this paper the probability density p on Q and the proba- 
bility current j on Q are permutation-invariant, i.e., p{x a \, . . . , x an ) = p(x\, . . . , x n ) and 
j*k(xai, x<m) = 3k{ x U ■ ■ ■ , x n) for any permutation cr of {1, . . . , n}, with j k G T X J1 
the component of j associated with the fc-th particle. 

We now define the relevant Hilbert space. The 1-particle Hilbert space associated 
with the time leaf E consists of square-integrable cross-sections of the spin bundle, 
if) : E — > U^gsS'a;, relative to the measure d 3 (x) dx. For this space we write 

3tf x = L 2 (E, U xeS S x , d 3 {x) dx) . (58) 

For any 1-particle Hilbert space Jt{, let T { (J^fi) and T h (M[) denote the fermionic and 
bosonic Fock space over M{, respectively; i.e., 

00 00 

r f (^i) = s_j*f n , r b (^) = s+jz** , (59) 

n=0 n=0 

where S_ and S + are the anti-symmetrizer and the symmetrizer, respectively. We write 
T x for either Tf or Tb- The relevant Hilbert spaces for us will be 

= r x (L 2 (E, U xe xS x , d 3 (x) dx)) (60) 

for E G T . A vector if) G can be regarded as a cross-section of the bundle U q S q , 
i.e., if) : <2 — > UqgQS'g with -0(g) G S^. In these terms, the inner product of can be 
expressed as 

(</#>= f dqd 3 (q)<f>*(q)if>(q). (61) 
J Q 
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Every density matrix p on can be expressed as a function p(q; r) of two ordered 
configuration variables. In fact, p(-; ■) is a cross-section of the bundle U qtr Hom(S r , S q ) 
over the base manifold Qx Q, where Hom(S r , S q ) denotes the space of C-linear mappings 
S r > S q . 

We are now ready to turn again to the quasi-Lindblad form as in (1181) . and to point 
out in which way ([I]) is of this form, provided we can leave aside the complications 
arising from the limit y — > x,y 5? and evaluate the density matrix directly on the 
singularity. Set A = (s,x) with s G {1,2,3,4} a spin index and x 6 S?t a point on the 
singularity; use the measure Pt = H z ® d 2 xw(x), where # denotes the counting measure 
on {1, 2, 3, 4}; that is, 

[ p t (dX)f(X) = V / d 2 xw(x)f(s,x). (62) 
J s=1 Js> t 

Let {bi(x), . . . ,6 4 (x)} be an orthonormal basis of S x consisting of eigenvectors of the 
positive definite Hermitian endomorphism 

ay(x) = lim c^x) d 4 (y)a^(y) (63) 



of S x with eigenvalues a\{x), . . . ,a 4 (x). Finally, let A t (s,x) be y/a s (x) times the an- 
nihilation operator on Fock space that annihilates a particle with position x and spin 

b s (x): 

(A t (s,x)^)(q) = y/^x)y/#Q + lK(x)ip(q,x). (64) 

Then 

(A t (s, x)pA* t (s, x)) (q; r) = ^ (#q + 1) (#r + 1) 6*(ar) p{q, x; r, x) a s (x) b s {x) , (65) 
and thus 

4 

J2{Ms,x)pA* t (s,x))(q;r) = V(#g + l)(#r + 1) tv Sx {p(q, x; r, x) a y (x)) , (66) 

3=1 

which brings ([1]) into the form (jl8p . 

4.2 Derivation From Equivariance 

We now give a derivation of §2§, and thus of the explicit evolution equation of the density 
matrix, from the framework of Bohmian mechanics, in particular from the requirements 
of equivariance and independence of disentangled systems, and from the assumption 
that the evolution of the density matrix is linear (i.e., of the form (pQ) with unspecified 
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1. We assume that particles move according to 

^=«i(t,Q(t)), (67) 

where v(t, q) is a time-dependent permutation-invariant vector field on configura- 
tion space <2(£t), and i runs through the 3#Q dimensions of Q(E t ) at Q. Here 
we use a coordinate system whose time component is given by the parameter t of 
the time foliation and that is otherwise arbitrary. (The vector field v depends on 
the choice of coordinates.) 

The motion according to fIBTj) continues until the configuration hits the boundary 
of configuration space, which means that one (or more) of the particles hits (or hit) 
the singularity; in this event that particle gets (or those particles get) annihilated, 
i.e., removed from the configuration, according to (fTUj) . 

2. We assume that 5?t = S t n=5^ is 2-dimensional and can be regarded as the boundary 
of S 4 . (In our Schwarzschild example, 



y t , = 5V n S = {f = const., x' = ±Vt' 2 - 2M] (6£ 



when t' > V2M; this is a disjoint union of two 2-spheres, and is indeed the 
boundary of £f, a disjoint union of two items of topology [0, oo) x S 2 .) 

Let vy(t,x) denote the speed at which the singular boundary J?t, regarded as 
a time- dependent 2-surface in coordinate 3-space, moves with increasing t in the 
surface- normal direction at the coordinate point representing x e =5^; since 5? 
is spacelike, this speed must be greater than the speed of light. This speed can 
be expressed in terms of the function T(£) that specifies the time coordinate at 
which a certain point £ in coordinate 3-space is reached by the singularity. (In our 
Schwarzschild example, if £ = <p) then T(£) = Vx' 2 + 2M.) Indeed, 

Vy(T(0,0 = ] ^77^, (69) 

|gradT(£)| 

where | • | is the Euclidean norm in coordinate 3-space and grad T the gradient of T 
in coordinate 3-space. To see this, note that a line in coordinate space-time starting 
at (T(£),£) on the singularity with direction (| grad T|, grad T/\ grad T\) will be 
tangent to the singularity. Note also that gradT(£) is orthogonal, in coordinate 
3-space, to ^r(£). (In our Schwarzschild example, gradT = x' ' / 'v 'x' 2 + 2M d x >, 
pointing in the radial direction, and Vj>(T = v 'x' 2 + 2M/\x'\.) 

It now follows that the continuity equation for a permutation-invariant distribution 
density p (relative to coordinate volume) of the configuration reads 

dp 3 * 9 
i=i 

+ (#9+1)/ d 2 xp(t,q,x) (v y (t,x) - v x , ± (t,q,x)) , (70) 



20 



where d 2 x is the surface area element in coordinate space (in our Schwarzschild 
example, d 2 x = dfidip), and v X) ± denotes the component of the velocity of the 
particle at x that is orthogonal, in coordinates, to 5?% and inward-pointing (i.e., 
away from the singularity). It can be expressed as 

gradT ( \ 

(^1,^2,^3)- (71) 



gradT I 



The first term on the right hand side of flTOj) . the negative spatial divergence 
of a probability current, represents the change in density due to the flow with 
velocity v(t,q), and the second term represents the gain in density due to jumps 
in configuration space to the configuration q from configurations containing one 
further particle at location x, where x lies on the singularity. To understand what 
the second term must be, keep q fixed, note that there are #g + 1 possibilities for 
the position of the variable x among the variables of q, consider configurations of 
the form (q, x) with x near the singularity, and note that the particle at x will be 
swallowed by the singularity within the next dt seconds if and only if its distance 
from the singularity, in coordinates, is less than (yy>(t,x) — v Xt ±(t, q, x))dt. Thus, 
the second term is the amount of p(t, q, x), for fixed q and arbitrary x, that flows 
into the singularity within the next dt seconds. 

3. We assume that the velocity vector field v(t, q) is given, as usual, as the quotient 
of a current vector field and a density function, 

^ = Mr (72) 

4. We anticipate that the current j and the density p are determined not by a wave 
function but by a density matrix p t , j(t, q) = j^(q) and p(t, q) = p^(q). Moreover, 
we assume that they are given in terms of pt by the formulas (fl2"j) and (114p . These 
formulas are the obvious extensions from wave functions to density matrices of the 
formulas (1251) and (I26I) usually utilized in Bohmian mechanics. In particular, for 
a pure state p = (H2D and dHJ reduce to ([25]) and fl26|) . 



5. We assume that the density matrix p t obeys, as in quantum dynamical semigroups, 
a linear evolution, which we write 

^ = {\p tt H t )+Cp t . (73) 

This is literally the same equation as (Q, but so far with £ unspecified. The 
operator H t is the Dirac Hamiltonian; we have separated this term knowing that 
it remains in the absence of the singularity J7, and thus know that £ = in the 
absence of a singularity. As a consequence of p^ evolves as follows: 

dv* ^ - 
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Here we have used the known continuity equation for the Dirac equation in a curved 
space-time, a key element of the proof of equivariance of Bohmian mechanics in a 
curved space-time [29, 30J. 

6. We want equivariance, i.e., pt = p Pt for all t if initially. For this, we need that 
the right hand sides of equations ( 1701 and ( 1721) coincide when p t = p pt is assumed; 
since the first terms (the divergence of the current) do coincide in that case, we 
only need that 

tr Sg (£p t )(g;g) = (# g + l)x 

x d 2 x lim trs q , y (pt(q,y]q,y)(d3(y)v,y(t,x) I - d^a 1 ^))) , (75) 

Jy t y^x,y<£y \ J 

where 

a±{y) = [g^ff ' a2{y) > ^ (76) 

is the component of a p {y) that is tangent to E f and orthogonal, in coordinates, to 
S?t at x. 

If (1751) is the case, we can argue that as soon as p t = p Pt for one t, then dp t /dt = 
dp pt /dt, and thus pt+dt — p Pt+dt , and so on into the future. Put more mathemat- 
ically, we take for granted that the continuity equation (ITUj) . as a PDE for p t (q), 
has unique solutions for every initial condition, observe that (t, q) i— > p pt {q) is a 
solution of (1701) by virtue of fl74l) and f!75l) . and conclude that if p agrees with p pt 
for one t then it must also agree at every later time. 

7. The obvious choice of £ that will make (!75|) true is 

(^)(g;r) = (#g + l) 1 / 2 (#r + l) 1 /2 x 

x/ rf 2 x lim tr s (p t (q,y;r,y)(d 3 (y)v^(t,x)I-d 4 (y)a ± (y))) . (77) 

The only differences between fl75l) and fl77|) are the trace over S q and that (1751) has 
q inserted for r. 

Indeed, (1771) is strongly suggested by the wish that disentangled systems should 
behave independently: Consider two disentangled systems, the first consisting of 
just one particle that is about to hit the singularity, the second consisting of several 
particles far away from the singularity. Then 

Pt{q, y\ r, y) = pj 1} (y; y) ® pj 2) (q; r) , (78) 

with the superscript indicating the system, and the contribution to pt with one 
particle less, Cpt(q;r), arising from the particle of the first system hitting the 
singularity, should be proportional to p\ (q; r). This is the case according to (J77J. 
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8. Eq. (J77I) agrees with (jSJ): Note that, in coordinates, 

1 — grad T 



with k = | gradT| 



w(x) 



VT 



K~ 



-2 ■ 



(79) 



JO) 



w(x) c^x) d 4 (y) ^(y) 



1 + K 2 K 
1 + K 2 K 



d±(y) a°(y) — d±(y) grad T • (a 1 , a 2 , a 3 " 



k 2 + 1 fd 3 (y) , j grad T , 2 



l + /t 2 



J - d 4 (y) 



(a , a , a ) 



This concludes our derivation. We have shown in particular that the distribution 
as in fll2p is equivariant with respect to the evolution of the Bohmian particles. 



In Schwarzschild space-time for t' > \JlM the evolution equation (JTJ) reads explicitly: 

dp t >(q;r) 



dtf 



tH t ,(q)p t ,(q;r) + ~H t >(r)p t >(q; 



+ AM^(#q + l)(#r + 1) V / dtf sin W , 



dtfsini? / rf<y9 lim r 3/2 (t' ,y') x 

y'\Vt' 2 -2M 



x tr s p t / (g, {oy\ cp); r, (try 7 , 0, y>)) 



Vt /2 - 2M 



/-F^a 1 ^'))! 
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with S y i the spin space at (t 1 , y', <£>) and 



F(t', yO^W) 



/0 1\ 

10 

10 

\1 0/ 



33) 



relative to the orthonormal basis of S y > associated with the orthonormal basis 

(F- 1 ^, F- 1 ^, r- 1 ^, (r sin^)- 1 ^) 



34) 



of T y u%. Note that all terms in ( 1821) of the form exp(r(t',y')/2M) could be dropped 
since r — > as y' — > \fV 2 ~^2M . Another remark concerns the term r 3//2 in (|82|) : This 
factor, which is essentially d 3 (y), tends to zero as y approaches the singularity; since 
the density in coordinates is conserved and thus cannot tend to zero when reaching the 
singularity, we must conclude that the density relative to the invariant measure, as given 
by trc pt(q, y; r, y), diverges at the singularity at the rate r~ 3//2 . 
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We close this subsection with a remark about the impossibility of a unitary evolution 
in the presence of a spacelike singularity. The process of the particle configuration has 
continuity equation (1701) . while a unitary evolution of a wave function, 

z^ = F V> + #W>, (85) 
at 

where Ho is the free Dirac Hamiltonian and H\ a putative further term, would imply 
for p^(q) = d 3 (t, q) iP*(q) if) t {q) that 

"if = " E + 2J ^ lm ^(q)(H^t)(q) . (86) 
i=i 

For equivariance we need this equation to coincide with ( 1701) . and using ( 1721) we find 
that we need that 

(#9 + 1) / d 2 xd 3 (x)4>t{q,x)if} t {q,x) (vy(t,x) - v x ,±{t,q,x)) = 



ft 



{Iv^rMiHiAM- (87) 



But the last equation cannot hold, no matter how we choose Hi, as the left hand side 
contains the factor ip^(q,x) but the right hand side has ipl(q) instead. In the case of 
timelike singularities, a similar problem can be solved by imposing a quasi-boundary 
condition on ip t that will ensure a relation between if>t{q) and ijjt{q,x); however, for a 
future spacelike singularity, there is no room for a boundary condition, since ipt{c[,x) is 
determined by the Dirac evolution. 



4.3 Quasi-Lindblad Equation for Reduced Density Matrix 

The quasi-Lindblad equation ([1]) with (T5]) also arises in another context, in which the 
singularity is replaced by a spacelike hypersurface 5?: then with (j2J) describes the 
time evolution of the reduced density matrix of what has not yet passed the hypersurface 
5? . So consider a space-time without singularities with a time foliation {£*} and therein 
a spacelike hypersurface 5? such that — E 4 D =5^ is always 2-dimensional. 
As a concrete simple example, readers may think of Minkowski space-time 

{JC, g) = (R\ diag(l, -1, -1, -1)) (88) 

and the time foliation defined by a Lorentz frame, so that S t are parallel spacelike hy- 
perplanes (where the time coordinate assumes the value t). Let 5? be another spacelike 
hyperplane, not parallel to the H t . We first describe the general definitions and then 
illustrate them using this example situation. 

Let if) : Usejc-E^ — > U£ e jc-U geS iv S q be a fermionic or bosonic iV-particle (Dirac) wave 
function evolving, for example, without interaction. Let J + (y) denote the future of 5? 
and J~{y) its past. For every £ e J 7 , set £± = £ n J ± (^). Let J? r± = {E 1 * 1 : E G J 7 }, 
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which is a foliation of J ± (y). It is a basic fact about (fermionic or bosonic) Fock spaces 
T x that 

r x (^ a © jt 16 ) = r x (^ la ) ® r x (jr 16 ) . (89) 

As a consequence, for disjoint subsets A, B of 3-space, A fl B = 0, 

T X (L 2 (A U B)) = T X (L 2 (A)) ® r x (L 2 (5)) . (90) 

Since S \ (S + U £~) = X fl 5? is a null set by assumption and thus not relevant to 
square- integrable functions, we have that, for every £ e JF, 

JT S = . (91) 

Now define 

p- =tr+|^)(^|, (92) 

where tr+ means the partial trace over ; p~ is a density operator on J^,- . 

We note that, as a consequence of the fact that ip lies in the A^-particle sector of 
Fock space, p~ is block-diagonal relative to the particle number sectors of Fock space, 
i.e., (4>\p~\x) = whenever lies in the n-particle sector of J^s- and \ in the m- 
particle sector with m ^ n. To see this, note that this is the case for |V's)(V's| by 
assumption, and the partial trace can be carried out using an orthonormal basis {bi} 
of that consists of basis vectors that are eigenvectors of particle number, so that 
(<i>\P~\x) = Ei(<f>®b i \fa)(fa\x®b i ). 

We can obtain a more explicit expression for p _ by writing it as a function p~(q; r); 
namely, with < n < N and g,r 6 (S _ ) n , 

p-(q;r)=r^J j dqd 3 (q) tr Ss (ip(q,q)ip*(r,q)) (93) 

(£+)iV-n 

with d 3 (q) the 3- volume factor as in fflBl . ip = and tr^, the partial trace over those 
spin indices belonging to particles in q. Here we use that 

S q ,g = S q ^Sg. (94) 

The binomial factor in (1931) arises from the re-ordering of variables in ip so that the 
N — n variables of q appear last. In case n = N, we take (E + ) Ar_n = (S + )° to be a 
one-element set and the 0-dimensional integration measure to give measure 1 to that 
one element, so that the integral equals ijj(q) if)*(r). In other words, the block pjj- of p _ 
in the N-particle sector of is just P\ifi)(ifi\P with P : L 2 {Y> N ) -> L 2 {{Yr) N ) the 
projection to the subspace in which all particles lie in E~. 

Let us formulate the time evolution of p~ . For an expression like dpj (q;r)/dt to 
make sense, we regard now q and r not as points on £ t but as their spatial coordinates; 
for simplicity, we will write E t for the image of H t in coordinate 3-space. We need 
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to differentiate fl93l) with respect to time, and thus to differentiate an integral with 
time-dependent domain, as in 



g(t,x) = / dyf(t,x,y) , (95) 

JBt 

where x is a variable in R™, y a variable in R m , and B t C R m a set with smooth boundary 
dB t moving in a smooth way. The rule we need can be regarded as a version of the 
fundamental theorem of calculus, which in its simplest form reads 

j t J*duf(u) = f(t) (96) 

while the form we need reads 

% = f d y~i + f d V v B(t,y) f(t,x,y) , (97) 
ot J Bt dt J 9Bt 

where VB{t,y) is the (signed) speed at which the surface dB t moves outward in the 
direction orthogonal to the surface at y. We thus obtain, with B t = (J2f) N ~ n , 

dpT(q;r) fN\ f dd 3 (t,q) 



dt \n J J dt 

+ f^j j dqd 3 {q) tr s .(-|^)(g, q) #(r, 0) + 
(E+F-» 

+ r) J dqd 3 (q) tr s . >fa(q, q) {(H t ^ t )*(r, q) + 

+ ( ^ J J dq d 3 (q) V B {t, q) tr 5 . ip t (q,q) ipt ( r > O) = 
a(s+) JV -™ 

= / 1 +/ 2 + J 3 + /4, (K 



where I\ + 12 + 1 3 correspond to the first term on the right hand side of fl97|) , and I4 to 
the second. 

To evaluate I±, note that dB t consists of N—n facets of the form (Sffl^ 7 ) x (E^)^" -1 
and permutations thereof; on the first facet, 

v B (t, q) = v^(t, q t ) , (99) 

where v y is the speed at which J?t — E t n =5^, regarded as a surface in coordinate space, 
moves in the normal direction. Exploiting the permutation symmetry of ip and 

<"-»>Q = (»+i)(„y (wo) 
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for n < N, we thus obtain 



l4 — (N — n)[ ) / dx d 3 (x) v y(t, x) x 



n 



x ^s x J dqd 3 (q) tr s . rj} t (q, x, q) $T(r, x, g) = 

( E +)N-„-l 

= (n + l) y dxd 3 (x) Vy(t,x) tr Sx p~(q,x;r,x) , (101) 

which is also true for n = N since then both sides vanish. 

The Hamiltonian in J 2 has two contributions, one acting on q and the other on q: 

H t = H t (q)+H t (q). (102) 

Correspondingly, we split 

I 2 = I 2 (q) + I 2 (q) . (103) 

The same calculation that leads to the continuity equation for the probability density 
and current of the many-particle Dirac equation also shows (when applied only to q) 
that the integrands of ii, h(q), and I 3 (q) together equal 



-^2diJi(q,r,q) := - ^ ^ d k ^(d 3 (q \ q k ) d 4 (q k ) tr s .(^ t (g, q) $ (r, q) a^{q k ))j , 

i=l k=l fi=l 

(104) 

where the ctj act only on the spin indices of particles belonging to q, not q or r. By the 
Ostrogradski-Gauss integral theorem, the integral of the divergence is the flux across 
the surface, so that, with n(q) the outward-pointing unit normal vector (in coordinate 
space) on the surface dB t at q, 

h + h{q) + h(q) = -() \ dqn(q) ■ J(q,r, q) = 

W J dB t 

= -{N-n)(^^j J dxd 4 (x)x 

x tr 5x J dq d 3 (q) tr s _ ip t (q, x, q) $\ (r, x, q) a x (x) = 

(E +)AT-„-l 

— —(n + i) J dxd^x) tr Sx f>t{q,x;r,x) a^{x) . (105) 



Here we used that on the first facet of dB t , i.e., on Sf t x (E^)^ - ™ -1 , the vector n(q) 
has 3(N — n) components of which only the first 3 can be nonzero, which form the unit 
normal vector on Sf t in coordinate 3-space. 
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In the term h(q), the Hamiltonian Ht(q) acts on a variable that is not integrated 
over, and thus can be exchanged with the integration, which leads to 

h{q) = -\H t {q)Pt (106) 

and likewise 

I 3 (q) = ip;H t (r). (107) 
Putting together ([98]), (jIUIj) . (jTQgj) . (1T05|) . ffTOBD . and flTOTj) . we obtain that 



dt 

with 



i[ft,H t ]+£p; (108) 



Cp t (q;r) = (#g + 1) 5 #9i#r ^ dx tr & (p t (g,x;r,x) (d 3 (x) x) I - d 4 (x) a ± (x)) ) . 

(109) 

If we drop the assumption that ip lies in the iV-particle sector of Fock space, we obtain 
instead of fl93l) that, for q,r G Q(£~), 

r(«r)= | dgd 8 (9) (^^^(^/^^^^-^^gl^^g)) (110) 

S(E+) 

and instead of (11091) that 
Cp; (?; r) = V(#?+l)(#r+l) x 

x / t t, s ,(«,, Iiril)Wl )v( 1 ,^-4W« i W)). (ni) 

This equation agrees with (I77|) and thus with (j2J), except that the need for a limit y — > x 
evaporates as the point x is now not singular. 

In our Minkowski example around (1881) . E$ is coordinatized as R 3 ; the Riemannian 
metric on £ f is the flat Euclidean metric on R 3 ; d^{x) = 1 and d±(x) = 1 everywhere; 
5* x = C 4 ; T(£) = T + i] ■ £ for the appropriate T G R and r\ G R 3 . We choose the 
Lorentz frame in such that To = and 77 = (0, 0, k) with constant < k < 1; such 
a choice is possible within those coordinate systems for which the time coordinate is 
constant on the time leaves. Now J/^, which is a 2-dimensional affine plane in R 3 , has 
the particularly simple form ^ = {£ G R 3 : £3 = t/n}; the gradient of T in R 3 is 
gradT(£) = 77; the speed at which 5^ t moves is Vy(t,^) = 1/| gradT(£)| = l/\r]\ = 1/k. 
Thus, 

£fh(q;r) = v/(#g + l)(#r + l) / d 2 x tr Sx (p7 (g, (x, £); r, (x, J)) (I - a 3A )) , (112) 
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where 



/O 







-1 



1 







"3,5. = 





1 



-1 





(113) 



0/ 



is the standard from of the third Dirac alpha matrix, acting on the spin index associated 
with the particle at (x,t/n). 

4.4 Absorbing Hypersurfaces 

We now present another derivation of our Bohm-type dynamics in the presence of a 
singularity, along the following lines. We first consider, instead of the singularity J?, a 
spacelike hypersurface J?, and set up a version of Bohmian mechanics in which every 
particle gets annihilated when hitting we call this model Bohmian mechanics with 
absorption at 5? . Then we let 5? approach a singularity 5? and argue that Bohmian 
mechanics with absorption at (the hypersurface) 5? converges to Bohmian mechanics 
with absorption at (the singularity) ,5^, i.e., to the theory described in Section [L41 

The definition of Bohmian mechanics with absorption at the hypersurface 5? is 
natural and straightforward. Let if) be the wave function as it would evolve without 
absorption, let p and J be the probability distribution and the probability current that 
if) defines on any E , let p and j be the appropriate marginals of p and j on £ _ , 
and use them to define Bohmian trajectories by means of the equation of motion pip . 
olQ/dt = j/p. Observe that p and j can be obtained directly, without reference to ip, 
from p~, which is defined as in the previous subsection as the partial trace of \if>)(ifj\ 
over the future of 5? and evolves according to the quasi-Lindblad equation fl 108[) with 
flllll) . In Bohmian mechanics with absorption at 5? , we do not mention if) but take a 
density matrix p t = to be fundamental and to evolve according to the quasi-Lindblad 
equation (11081) with (JUT)) , and we define p, j, and the law of motion in terms of p t . We 
now give the defining equations. 

The probability density p(q) on configuration space <2(£~) for £ G T is defined in 
the following way. Let p be the probability density that would arise from if) without 
absorption at p is defined on Q(S). Using the identification 



if A fl B = (consisting of re-ordering a configuration q in such a way that the particles 
in A are listed first, q = (qnA,qnB)), and ignoring £ \ (S + U £~) because it has lower 
dimension, the function p can be written as p(q, q) with q G Q(S~) and q G Q(S + ). We 
have that 



Since we assume that particles vanish when hitting the hypersurface 5^ while other 
particles are not affected, we are led to 



Q(AUB) = Q(A) x Q(B) 



(114) 



p{q, q) = d 3 {q, q) ip*(q, q) ip(q, q) . 



(115) 




(116) 



Q(E+) 
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This can be re-expressed in terms of p , see f l93l) . as 

p(q) = d 3 (q) tr Sq p~ (q; q) , (117) 

which is the density in configuration space associated with p~~ in the natural way, parallel 
to pi). 

Likewise, the probability current vector field j on configuration space is defined as 
follows. Let j be the probability current vector field obtained from if) and defined on 
Q(E): for q G (£~) n and q G (£+) m , 

fk(q, q) = dz(q \ q h ) d 3 (q) d 4 (q k ) ^*(g, q) a»(q k ) 1>(q, q) (118) 

with k — 1, . . . , n. We define j to be the marginal current on Q(E~), or, in other words, 
the average current given q; explicitly, for q G (£~) n , 

#(?) = J dqf k (q,q). (119) 

s(s+) 

As a consequence, 

= 4(g \ g*) 4(g fc ) tr 5 , (p- (g; g) a M (g fe )) • (120) 

Now take (11081) with fillip to be the fundamental evolution law of a density matrix 
p t = fa . For these equations to define an evolution law, we need that is spacelike0 
which we always assumed, and the following locality property of the one-particle Dirac 
equation^] Suppose is a one-particle wave function; in order to predict <f)(x) for some 
space-time point x, one needs only the initial data in the past of x, i.e., 0|s o nj-(x) if 
Eo is the surface on which the initial data are specified. In particular, if x G J~{y) 
then S fl (x) C (S") . As a consequence, the time evolution equation (11081) can be 
solved on J~(J?) without ever computing ip at any configuration containing any point 
from J + (^), and without knowing the electromagnetic field or the metric at any point 
of J + (y); indeed, the evolution equation (11081) can be solved using initial data pg on a 
surface E that already intersects 5? . 

Having obtained the density matrix p t = p^ , let it define p and j as above, and let 
the particles move according to (|lip . whenever none of the particles hits the absorbing 
surface y. As soon as a particle hits 5? y that particle gets annihilated and removed 
from the configuration, as in (fTOj) with replaced by 5? . 

Let us turn again, in more detail, to the difference between ignoring and absorbing 
particles behind 5? . While the density p and the current j on S~ are merely the 
marginals of the density p and the current j on S, the Bohm-type trajectory Q t obtained 

5 In case 5? is not spacelike, as long as it divides ^ into two connected components ^# ± , i.e., 
J(\ 3* = u J?-, then, for Y.f = E t n ^# ± , the derivation of ([TU5|) with (fT09|) remains valid, but 
does not uniquely determine the evolution of from an initial datum p~[ . 

6 Most other relativistic wave equations share this locality property, in particular the Maxwell and 
Weyl equations, but not the first-order Klein-Gordon equation idtip = yrn^— Atp. 
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from p and j by ffTTj) are usually very different from the trajectory Q t obtained from 
p and J (and thus from %p): Whereas the velocity of particle 1 in the configuration Q t 
may depend on the position of particle 2 that has already crossed 5? , the velocity of 
particle 1 in Q t does not so depend; after all, particle 2 does not exist any more after 
hitting 5? ', and its position is therefore not defined. Instead, the velocity of particle 1 
in Q t is equal, when expressed in terms of ip, to the average (of the velocity of particle 
1 in Q t ) over all positions that particle 2 might assume behind 5? . Hence, while on one 
mathematical level — the level of wave functions, density matrices, probability densities 
and currents — it may seem like our construction merely involves ignoring the particles 
behind on another mathematical level — the level of the trajectories — the absorption 
of particles at has an effect on the other particles that mere ignoring would not have. 

The evolution of the Bohmian particles just defined is equivariant, i.e., if the config- 
uration Q t on E t is random with distribution density p^ then for every s > t, Q s has 
distribution density p ps . This can be shown with the same argument as used for (170|) . 
flU}, and CED. 

Now we assume the existence of a spacelike singularity 5? in the future of 5? . Then 
the evolution of ip, from which we obtained that of p~ by a partial trace in Section H~3"| 
is not defined any more, but the evolution (11081) of p~ is still defined because of the 
locality property mentioned above. As we let the spacelike hypersurface approach 
=y, (I108p with (111 II) formally converges to ([1]) with (jSJ), while the laws for the Bohmian 
configuration remain unchanged. We thus obtain the Bohm-type evolution as in (fTT)- ffTUl) 
in a different way, as a limit of the evolution with absorption at a hypersurface. 

4.5 Abstract Mathematical Structure of the Time Evolution 

The time evolution we are considering, as summarized by flTJ), does not fit into the 
usual categories of quantum mechanical time evolution, as it corresponds neither to a 
unitary one-parameter group U t on a fixed Hilbert space nor to a quantum dynamical 
semigroup (as would arise from a Lindblad equation) on M' . Rather, with every time 
leaf Et there is associated a Hilbert space J^, and the time evolution from E s to E$ 
corresponds to a superoperator 

S\ : TRCL(Jf s ) -> TRCL(J^t) , (121) 

where TRCL(Jif) means the trace class of J$? (roughly, the set of operators Jif — > 
with finite trace). 

Let me elaborate a bit on the fact that Jif s and M[ are not the same space: It is 
always the case in relativity that the time evolution takes place from one hypersurface 
E s to another E t , and therefore that containing functions of configurations in E s , is 
different from containing functions of configurations in E t . On the other hand, one 
may often seek a way of identifying E s and E 4 , for example by means of coordinates. In 
the presence of a spacelike singularity, however, coordinates x M for which x° is timelike 
and spacelike may have the feature, like the Kruskal coordinates in the 

Schwarzschild space-time, that E s and E 4 correspond to different subsets of coordinate 
3-space, so that the coordinates do not provide a diffeomorphism E s — > E t . 
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So we are forced, more than ever, to regard and 3^ t as different spaces. As a 
consequence, the time evolution superoperator S\ must have two indices, indicating the 
initial time s and the final time t. And as a consequence of that, they cannot form a 
semigroup, since a semigroup {g t } requires that g t g u = gt+u- Instead, the appropriate 
notion of time evolution for our purpose is that of cocycle, a notion taken from the 
theory of random dynamical systems [27J and designed for describing the evolution in 
the presence of noise or time-dependent external fields. 

We define a quantum dynamical cocycle to be a 2-parameter family of mappings 



for < s < t such that 

• Si is C-linear, completely positive, and trace-preserving; 

• S?Sl = <S S " for < s < t < u; 

• Sl{p) = p for every p G TRCL(J%); 

It is plausible that equations (jTJ and (TSJ) define a quantum dynamical cocycle 5*: To 
begin with, (CQ) is linear in p t . To argue that Si is completely positive, we note that the 
composition of completely positive superoperators is completely positive; regarding S f s 
as the composition of many 

S u+du = i[-,H t ]dt + £ dt, we need that both the Hamil- 
tonian evolution and £ are completely positive, which is plausible. The conservation of 
tr p t corresponds to the conservation of total probability, which is exactly what £ was 
designed for. 

Another requirement that one may wish to add to the definition is that Si depend 
continuously on s and t, but for this one needs a topology on the bundle U t J4ft, which 
we have not defined yet. More importantly, when requiring continuity in s and t then 
the definition fails to cover the cases in which £ t H 5? is 3-dimensional, see ( 11291 ). 

Note that a quantum dynamical semigroup forms a special case of a quantum dy- 
namical cocycle, in which all Hilbert spaces Mi are identified with one fixed space ffi 
and S\ depends only on the time difference t — s. 

Another special case arises for a space-time without spacelike singularities, in which 
the quantum-mechanical time evolution is unitary. Then <S*(p) = UgpUf, where the 
£7* : J4? s — > J4?t are unitary isomorphisms implementing the unitary evolution of the 
wave function according to if) t = U^s] they satisfy Ufa = ip an d U^JJ\ = for all 
s,t,u G KL In this case we call (5^)o< s <t a unitary cocycle. As an example of a unitary 
time evolution that does not correspond to a unitary 1-parameter group, but can be 
represented as a unitary cocycle, consider a time- dependent Hamiltonian H t on a fixed 
Hilbert space, as in non-relativistic quantum mechanics with time-dependent external 
fields; then (for bounded H t ) 



S l s : TRCL(Jf s ) -> TRCL(J%) 



(122) 




(123) 
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with Texp the time-ordered exponential; in other words, Ul is the solution of the initial- 
value problem 

dU* - - 

-^ = -lH t Ul U: = I. (124) 

In order to define a Hamiltonian operator, or infinitesimal generator, for more general 
unitary cocycles for which the J% are not identified, we need a differentiable structure 
on the bundle U t J%t of Hilbert spaces. We begin with outlining an infinite-dimensional 
analog of the concept of a vector bundle. 

By a Hilbert bundle over the interval (ti,t 2 ) CRwe mean a family of Hilbert spaces 
J#t, t G (£1,^2), together with an equivalence class of bundle maps; a bundle map or 
trivialization is a family V t : J$ft ~ * ^ of unitary isomorphisms to some fixed Hilbert 
space J^; two bundle maps V t , V[ : — > M" are equivalent if and only if for every 
■0 G and one (and thus every) unitary isomorphism W : J^" — > the curve 
i i-> WVjVf 1 ^ in is We simply write U t J% to denote the bundle (with the 

understanding that J#l H = for t 7^ s), and V" : U t J£t — > (ti, t 2 ) x to denote the 
bundle map consisting of the Vt : M[ — > J^ 7 . 

The definition of a Hilbert bundle over manifolds other than intervals would in- 
volve local bundle maps (or local trivializations) and requirements on the transition 
functions between two local bundle maps concerning smoothness and consistency. The 
1-dimensional case, all we need here, is particularly simple. In our setting, the base man- 
ifold of the Hilbert bundle is really the time foliation J 7 , regarded as a 1-dimensional 
manifold containing one point representing each S t ; since we parameterize the time 
leaves E t by real numbers t anyhow, we can afford to regard the base manifold as simply 
an interval (ti,t 2 ) of t-values. 

Now we want to define an extension to Hilbert bundles of the concept of a connection 
on a vector bundle. A connection essentially amounts to an identification between the 
fiber spaces over infinitesimally neighboring base points, leading to a path-dependent 
identification between the fiber space over any two base points, known as parallel trans- 
port. Relative to a choice of bases in the fiber spaces, the connection can be expressed in 
terms of connection coefficients T a ba , the best-known example of which are the Christof- 
fel symbols, the coefficients of the connection naturally associated with a Lorentzian (or 
Riemannian) metric. Of the three indices a, b, a, the last one refers to the base manifold 
(for the Christoffel symbols, space-time; for us, the time axis as represented by J 7 ), 
while a and b refer to the fiber space (for Christoffel symbols, again the tangent space to 
space-time; for us, 3%t). Since for us, the base manifold is one-dimensional, we can drop 
the index a, taking d/dt as the basis vector along the base manifold; what remains is 
T a b , an operator on J#t, which we write as —iH t /h; if parallel transport respects inner 
products, then r a b must be skew-adjoint, or H t self-adjoint. 

We can thus define a connection on a Hilbert bundle U t Mt over (ti,t 2 ) to be given 
by an equivalence class of triples (Jf? , V, {H^)te{t\M))i wnere ^ is a Hilbert space, 
V : UtJ^t — > (ti, t 2 ) x is a bundle map, and (H t ) te ^ 1)t2 ^ is a 1-parameter family of 
operators H t on Jtif. If every H t is self-adjoint then we call the connection unitary. Two 
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triples (J^ 7 , V, (Ht)) and (J&", V', (H' t )) are considered equivalent if and only if 



UvlY^vX- 1 = H t - ihjiy^vl)- 1 )^ 



(125) 



where the derivative is understood as a strong limit. To understand this definition, think 
of the connection as defining a parallel transport operator [/* : J^f s — ► by means of 



Ul = V^Texp^-ih" 1 J H u dujv a . 



(126) 



Since U\ should be independent of the bundle map V, H t has to transform in the 
appropriate way, and a short calculation shows that (11251) is the transformation law. 

As an example of a Hilbert bundle and a connection, consider a single Dirac particle 
on a singularity- free space-time j% with global coordinates (x°, . . . , x 3 ) : ^ — > M 4 , 
forming a diffeomorphism, such that x° is timelike and spacelike. Then 

E t = {a; = t} is a spacelike hypersurface, and ^ = (x 1 ,^ 2 ,^ 3 ) : E 4 — > M. 3 is a 
diffeomorphism. The Hilbert space M[ consists of measurable cross-sections of the spin 
bundle U x£ s t S x that are square-integrable relative to the Riemannian 3- volume measure 
on E t and relative to the inner product on S x associated with the unit normal vector on 

at x. Let J$f = L 2 (R 3 , C 4 ) with respect to Lebesgue measure on M 3 and the standard 
inner product on C 4 . Note that the Riemannian 3-volume measure has non-constant 
density ds(x) = a/— det g^(x) relative to the Lebesgue measure on coordinate space 
IR 3 , so we need to compensate for that. At every x G E t choose an orthonormal basis b x 
of S x that depends smoothly on x. Obtain b x from b x by scaling each basis vector by 



a/g? 3 (x). Then define V t by 



(V^) s (q) = b* x Mx) 



(127) 



where q G M. 3 is a point in coordinate space, s G {1,2,3,4} a spin index, b x>s the s-th 
element of b x , and (f>*tp the inner product in S x . To see that Vt is unitary, note that 



dq^2<f)*(x)b x>s b* XjS i/;(x) 



x =<p t (?) 



dq d 3 (x) 4>*(x)ip(x) 



and that V% is clearly surjective. The Dirac equation defines a unitary time evolution tP a : 
E s — > E t , corresponding to a unitary connection on U t J^, expressed in the coordinates 
x M as a time-dependent self-adjoint Hamiltonian i^. 

In the presence of a future singularity, the J% still form a Hilbert bundle. For 
example, in the Schwarzschild space-time with the t'-foliation for t' > and J% = 
r x (L 2 (E t /, *S f |s t /)) ; where T x means either the bosonic or the fermionic Fock space, we 
can set J$? = T X (L 2 (R x § 2 , C 4 )). For < f < 2M, we use the coordinates as before 
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to identify with R x § 2 . For i' > 2M, £f has two connected components, E^, with 
x' > \Jt 12 — 2M and £j, ' with x' < — \Jt' 2 — 2M ; by replacing the x' coordinate on E^ 
with x" = x' =F yt/ 2 — 2M, we identify E# with (R \ {0}) x § 2 ; the one point missing 
from R does not affect the L 2 space. Then the bundle map V can be defined as before. 

In the presence of a future singularity, the Dirac equation fails to define a unitary 
evolution, but it still defines a time evolution for a wave function <p t e J#h t , one for which 
\\<f>t\\ decreases with increasing t. For example, in the setting of Section involving 
instead of a singularity an absorbing spacelike hypersurface y , (f>t is just the restriction 
of ip t (the wave function in the absence of the absorbing hypersurface) to Q(E^), i.e., <f> t 
is ip t evaluated only at configurations for which all particles are located in Ejr, i.e., in the 
past of the hypersurface correspondingly, ||0t|| 2 is the probability that no particle 
has hit 5? up to time t. In other words, in the presence of a future singularity the Dirac 
equation still defines a connection on U t J%t, the Dirac connection, but it is not unitary; 
correspondingly, the Dirac Hamiltonian is not self-adjoint. 

In terms of the Dirac connection, the evolution (TjQ) can be expressed as 

Vp t = Cpt, (128) 

where V is the covariant derivative operator associated with the Dirac connection (and 
applied along the vector field d/dt on the time axis), and p t is now a cross-section of 
the Banach space bundle UtTRCL(J#t) arising from the Hilbert bundle UtJtft. 



4.6 Remarks 

A Chunk of Singularity. In equation (pQ) we assumed that E t n 5? is 2-dimensional; if 
it is 3-dimensional for a particular t then p t should have a discontinuity as a function 
of t at t , according to 

Pt Q +{q; r) = p to -{q; r)+ J dq tr s . p to -{q, q\ r, q) (129) 

S(s* n,r) 

for q,re Q(E to \^). 

Non-Smooth Singularities. Penrose [25] has suggested that the future singularities 
arising from gravitational collapse may be rather irregular, and this further suggests 
that they may in fact be non-smooth. That is, after a conformal transformation the 
singularity may correspond to an achronal surface that is not smooth but merely con- 
tinuous. In our discussion so far we assumed smoothness, but it seems plausible that 
this assumption is not needed. 

The reason for believing this is that smoothness is not needed when considering an 
absorbing spacelike hypersurface y instead of a singularity 5? . Suppose 5? is merely 
continuous. If 5? fl E 4 is a null set in every E t , then Ef and pf are still well defined, 
since their definition did not involve differentiation. That pj is well defined is all we 
needed to show. 
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To be sure, p~j~(q,r) may fail to be differentiable with respect to t as a consequence 
of the lack of smoothness in the t-dependence of Y,f; however, pj may also fail to be 
differentiable with respect to t if any of its eigenfunctions is not contained in the domain 
of H t . After all, the Schrodinger equation 

ih^- = H t il) t (130) 

holds literally only for ip t in the domain of H t ; if ip to lies outside the domain of H to then 
t i— > ip t is not differentiable at t . To sum up, even though fll08j) cannot be expected to 
hold literally, the evolution of p^ should exist for non-smooth 5? . 

Objections. Banks, Susskind, and Peskin [2J have argued, in response to Hawking's 
[T8] proposal that the fundamental time evolution might transform pure to mixed states, 
that such an evolution would have to either allow superluminal signalling or violate 
the conservation of energy-momentum. Indeed, in our model energy-momentum is not 
conserved, as the energy- momentum of a particle hitting the singularity gets lost; it 
should be accounted for by a suitable change in the space-time geometry, but our model 
does not do that. On the other hand, our model does not allow superluminal signalling. 
This is most directly conveyed by the parallel between the evolution equations of the 
density matrix p t in the presence of a future singularity and the density matrix p^ = 
tr + \ip)(ip\ as m f|92|) obtained by tracing out the future of a spacelike hypersurface in a 
non-singular space-time (see Section and interpreting the operation of tracing out 
as merely ignoring part of the information encoded in if). It is known that the unitary 
evolution of ip does not allow superluminal signalling, and rather clear that ignoring 
some information cannot create a possibility of superluminal signalling. Still, it would 
be desirable to have a carefully formulated no-signalling proof. 

Maudlin [23J has argued that the pure-to-mixed evolution be an artifact of consid- 
ering the wrong spacelike hypersurfaces. If E s , he argued, is a Cauchy surface and E t is 
not, for example if they are hypersurfaces of constant t! in the Schwarzschild space-time 
with -y/2M < s < \p2M and t > \p2M, then it is no wonder that the evolution from 
£ s to S t is pure-to-mixed: after all, if S 4 is not a Cauchy surface then it is not adequate 
for describing initial data of the evolution. Thus, the pure-to-mixed evolution does not 
mean that anything is unusual about the evolution but that some hypersurfaces are 
inadequate. For example, the hypersurface {t = 0, x > 0} in Minkowski space-time is 
spacelike but not Cauchy — it is too small — and so the quantum state associated with 
it is the density matrix arising from the wave function on {t = 0} by tracing out the 
degrees of freedom associated with x < 0. If one sticks to Cauchy surfaces, then the 
evolution remains unitary, and black hole evaporation never occurs. However, this ar- 
gument does not work in Bohmian mechanics as it ignores the role of the time foliation. 
As mentioned at the end of Section 11.31 the time leaves S e T may not be Cauchy 
surfaces but instead border on the singularity. 

Black Hole Evaporation. There is no consensus in the literature about whether 
information is lost during black hole evaporation, i.e., whether unitarity is violated. 
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While in our model unitarity is indeed violated, this does not support conclusions about 
black hole evaporation, for two reasons: First, we assumed that the gravitational field 
can be described by a classical Lorentzian geometry with a spacelike singularity. This 
assumption we might be violated, as it may be necessary to apply a quantum gravity 
theory, and it can be questioned whether then any singularity will actually arise. Second, 
our model ignores any back reaction of the particles on the space-time geometry. It is 
clear that a black hole will grow in mass when swallowing particles, so the quantum 
state of the gravitational field should be affected by the infalling particle and might 
store the information lost from the quantum state of the matter, leading to unitarity of 
the full evolution of both matter and the gravitational field. 



5 Past Spacelike Singularities 

Our method of studying past singularities is to postulate reversibility of the fundamental 
laws of the theory. In this way, the laws we already have for future singularities determine 
the laws for past singularities. 

The obvious fact about past singularities is that no future-directed timelike curve 
can end there. Thus, a past singularity can emit but not absorb particles. Examples 
of past singularities include the set in Schwarzschild space-time (see Section [3]) 
corresponding to a white hole and the big bang singularity in Friedmann-Robertson- 
Walker space-times [19|. 



5.1 Evolution of the Density Matrix 

The equation (pQ) for evolving a density matrix to the future in the presence of a future 
singularity can be used, when time-reversed, for evolving a density matrix to the past 
in the presence of a past singularity. The time-reversed form reads 

0f>l -Up t ,H t ) + Cp t (131) 



d(-t) 

with C defined as in Compared to the equation ([T]) for a future singularity, the term 
Cp t has the opposite sign. If we choose an initial wave function on a Cauchy surface 
(such as {£' = 0} in Schwarzschild space-time), the density matrix p t is defined on every 
time leaf Ej. 

However, the time evolution towards the future is not well defined in the presence of 
a past singularity (like the time evolution towards the past in the presence of a future 
singularity): We have to invert a pure-to-mixed evolution, and since this evolution is 
many-to-one, its inverse is not unique. The problem is analogous to that of recovering 
a vector ip e M{ ® from its reduced density matrix p 1 = tr 2 \ip)(ip\. Thus, the 
evolution of p t towards the future is not uniquely determined. It is not governed by a 
stochastic law, either, but such a law could be added as follows. A theory could provide 
a probability distribution p^idipT.) for the wave function ^ s on a Cauchy hypersurface 
E, preferably in a way that does not single out any particular E. (For examples of 
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probability distributions over wave functions, see [IS].) Then any initial datum p to on 
a time leaf E to bordering on the past singularity (and thus in the past of E) defines a 
conditional distribution p-£(dips\Pt ), concentrated on the set of those ips's which, when 
evolved backward to E to , lead to p to , and thus also defines a stochastic process (pt)t>t - 

5.2 Evolution of the Configuration 

Let us ignore the problem of finding the density matrices pt in the presence of a past 
singularity — let us suppose we are given all density matrices pt — and focus on how to 
define the evolution of the configuration Q t . Sticking to reversibility, we obtain, from 
the evolution we know in the presence of future singularities, that the particles move 
according to ([7]), the Bohm-type law of motion using a fundamental density matrix, 
with new particles created at the singularity in a stochastic way, given by (TTOT) . 

Indeed, think of the evolution in the presence of future singularities as a stochastic 
process, i.e., as a measure P on path space. This process is in fact deterministic, that 
is, the initial configuration Q to is random with distribution p pt o, and the path is a 
function of p to and Qu. The time reversal mapping T maps every path to its time 
reverse, and P to T^Pjj Our claim is that T*P corresponds to a Markovian stochastic 
process with particle creation at rate (FIT?]) . Being defined by a probability distribution 
on path space, it obviously is a stochastic process. The Markov property follows from 
the determinism in the opposite time direction: If the past path is a function of the 
present configuration (for fixed p to ), then conditional probabilities of future events given 
the past path equal those given the present configuration. Since between two jumps 
the trajectory in configuration space is deterministic in both time directions, the only 
randomness concerns when to jump and where to jump. The only possible jumps are, 
up to permutation of the configuration, 

Q(t-)^Q(t+) = (Q(t-),X) (132) 

with X G corresponding to the creation of a new particle at the past singularity. 
(The creation of two or more particles at the same time has probability density zero.) 
To determine the rate of such a jump, note that in the other time direction, 

F(Q(t) G dq x E t , Q(t + dt) G dq, end point G d 2 xj = 

p(t, q, x)dq d 2 x (vy(t, x) — v Xi ±(t, q, x))dt , (133) 

and therefore 



TJPf Q(t) G dq, Q(t + dt) G dq x E t , creation point G d 2 x 




pit, q, x)dq d 2 x \Vy(t, x) — v X) ±(t, q, x))dt , (134) 

7 In the general relativistic context, time reversal is essentially a trivial operation, since the time 
reverse of a space-time, decorated with world lines, is isometric to the original, and thus physically 
equivalent. Still, a theory may fail to be reversible if it assumes, as we did, a time orientation. 
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which implies 



a t {d 2 x\q) dt = T*P ^creation within time dt and location d 2 x Q(t) E dq 
= (#g + l ) P ^^ d 2 x {vy(t, x) - v X)± (t, q, x))dt 

lim „ tT s q , y Pt(q, V, q, y) (d 3 (y)v^(t, x)I- d A (y)a x {y)) 

= [#q + 1) — r d xdt , 

tis q pt{q;q) 

(135) 

which agrees with (fT9l) . 

This evolution of the configuration, based on the combination of (JTj) and ( fl9l) . is 
equivariant, i.e., if Q(t) is random with distribution density p pt defined in (TT2l) then also 
for every s > t, Q(s) has distribution density p Ps . This follows from the fact that the 
time-reversed process is equivariant. 

5.3 Comparison with Bell- Type Quantum Field Theory 

Stochastic jumps of the configuration also occur in an extension of Bohmian mechanics 
to quantum field theory (QFT), known as Bell-type QFT [HJdHII]. Let us compare the 
jumps in the two theories. 

In Bell-type QFT, the jumps usually correspond to creation or annihilation of parti- 
cles, while a past singularity can create but not annihilate particles (assuming that the 
world lines are causal). Also in Bell- type QFT, the configuration process is Markovian, 
with the jump rate specified by a law |TT] that reads (after replacing the wave function 
in Fock space by a fundamental density matrix on Fock space) 

= (136) 

with Hj the interaction Hamiltonian, x + = maxjx, 0} the positive part of x G R, and 
P(dq) a positive-operator-valued measure (POVM) on configuration space serving as 
the configuration observable, usually 

P(dq) = \q)(q\®I Sq d 3 (q)dq. (137) 

Both (H2J and (ITHBjl are of the form 

tr(p t P(dq)) 

but with different operator- valued measures R(-) on Q x Q: In Bell- type QFT, 

R(dq' x dq) = -^P{dq')H I P{dq) , (139) 
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whereas in our jump rate (fl9~j) . 
R(dq' x dq) = (#g + 1) / P(dq x d 2 x) w(x) lim c M (x) d 4 (y) 6 ((Z)a:) (cV) (140) 

with 5 q the Dirac delta measure centered at q. Of course, since f)140p is positive-operator- 
valued, in this case the operations in (11381) of taking the real part and the positive part 
are trivial. 

The common structure (I138P is owed to the fact that both rate formulas can be 
obtained starting from the appropriate formula for the probability current 

J t {dq', dq) = Re tr(p t R{dq x dq)) (141) 

between volume elements dq' and dq in configuration space. The minimal jump rate 
<Jt(dq'\q) compatible with this current is (I138p . The form of the operators R(dq' x dq) 
is (not uniquely determined but) suggested by the probability balance equation ( 1741 . 
respectively the probability balance equation of QFT, which we need to agree with the 
probability balance equation of a jump process, 

i_1 <?'es 

where j is the current due to continuous motion and J the current due to jumps. 
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